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FOREWORD 


This  is  the  first  annual  technical  report  prepared  and  submitted  by 
Systems  and  Applied  Sciences  Corporation  under  Contract  F19628-80-C-0058 . 
The  report  discusses  the  development  of  analytical  techniques,  and  their 
implementation  through  computer  routines,  for  the  reduction  and  analysis 
of  digitized  data  transmitted  from  sensors  flown  on  research  balloons 
and  on  research  rockets  in  the  50  -  500  km  altitude  range.  The  data 
transmission  systems  employed  permit  the  use  of  a  standard  procedure, 
regardless  of  the  experiment  flown.  This  report  details  that  standard 
procedure  and  gives  an  example  of  one  experiment. 

Work  was  performed  for  the  Analysis  and  Simulation  Section  (SUWA) , 
Computation  Branch,  Office  of  Research  Services,  Air  Force  Geophysics 
Laboratory,  Hanscom  AFB,  MA  01731. 

SASC  technical  personnel  employed  on  the  contract  were: 

Paul  F.  Hilton,  M.S.,  Principal  Investigator 
Charles  F.  Ivaldi,  Jr.,  B.S. ,  Programmer 
Barry  A.  Mareiro,  Jr.,  Junior  Programmer 
Stephen  C.  Phillips,  Junior  Programmer. 

Program  management  was  provided  by  Alan  M.  Gerlach,  Ph.D. 
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A.  INTRODUCTION 


The  Analysis  and  Simulation  Section,  Computation  Branch,  Office  of 
Research  Services,  Air  Force  Geophysics  Laboratory,  i rovidos  a  wide  rancje 
technical  support  to  in-house  scientists  in  such  areas  as  staiist'Cal 
and  numerical  analysis,  computer  simulation,  and  reduction  and  analysis 
of  atmospheric  and  near-earth  space  environmental  data.  The  maic  r  e‘  ‘"orts 
are  in  those  areas  of  qeophyrical  research  which  utilize  data  tians- 
nitted  from  sensors  flown  on  research  ballcons  and  rochets.  The  respon¬ 
sibility  of  Systems  and  Applied  Sciences  Corporation  (SASC)  is  to  provide 
support,  not  available  from  in-house  sources,  in  the  reduction  and  analy¬ 
sis  of  rocket-  and  balloon-borne  sensor  data. 

Reduction  of  scientific  data  requires  a  wide  variety  of  analytical 
methodologies,  such  as  numerical  techniques,  including  averaging,  smooth¬ 
ing,  and  filtering  of  data;  ana  the  use  of  {polynomial  and  Fourier  approx- 
mations  to  model  the  data.  Generally,  the  analysis  is  carried  out  on  the 
AFGL  CDC  6600  computer  system.  Software  development  as  reauireu  fc'r  this 
contractual  effort  is  the  responsibility  of  SASC  technical  personnel; 
occasionally,  existing  routines  are  obtained  from  the  AFGL  oomputer  sys¬ 
tem  library. 

The  refined  scientific  data  are  correlated  with  trajectory  and 
attitude  information  to  provide  space-time  loci  for  the  data  points, 
and  are  displayed  according  to  the  specifications  of  the  experimenter. 
Display  methods  include  2-dijneu3ional,  3-dimensional,  and  contovn.  vilots, 
drawn  with  a  pen-and-ink  drum  {iotter,  microfiche,  or  Tektronix. 

The  various  research  areas  which  SASC  supports  may  be  grouped  into 
f'wo  major  categories:  particles,  plasmas,  and  fields  in  the  ionosi)he''f.  and 
magnetosphere;  and  optical  radiations  (infrared,  visible,  ultraviolet) 
from  auroral  and  other  natural  atmospheric  sources,  and  from  man-made 
sources.  Instruments  flown  characteristically  include  radiometers,  pho¬ 
tometers,  circular  variable  spectrometers,  electrostatic  analyzers, 
paiticle  counters,  electron  spectrometers,  energy  deixjsition  scintilla¬ 
tors,  and  plasma  frequency  probes. 

By  means  of  these  kinds  of  instruments  and  other  sensors,  experi¬ 
ment.''  are  carried  out  to  meet  a  wide  range  of  scientific  objectives.  For 
example,  radiometors  are  designed  for  measurinq  selected  auroral  emissions 


ill  ci'.ti  ii'ifrai.ed  reqion  thu  si-o.;i.ruii!,  v.’h:i.'  ;4i')t.!tT.L;ti_i  l;  i.  ci;  ;  o; 

moasurement  of  visible  and  ultraviolot  ein.ssions.  i'artioii.-  {'oanter s  ;.ro- 
vide  quantity  and  low  resolution  energy  mt-asure.-nent.-,  o:  enerqetic  elec¬ 
trons  within  a  predetermined  field  of  view.  h'ith  the  aid  of  plasma  fre¬ 
quency  probes,  electron  densities  and  te.mperatures  can  be  determined  ty 
resonant  frequency  techniques. 

Although  there  are  many  areas  of  research,  there,  i  r,  gent  roll  /  or.'.- 
data  reduction  and  analysis  methodology,  which  consists  of  fear  .iiscreto- 
stepa:  (1)  vehicle  profiles,  (2)  telemetry  data  with  standard  processing, 
(3)  merging  of  trajectory  and  attitude-  data  tapes,  and  '4)  standard 
reduction  and  display.  The  next  section  will  describe  the  standard¬ 
ized  data  processing  system  and  the  interactions  between  SLiWA  and  .SASC. 

B.  PROCESSING  SYSTEM  FOR  ROCKET-  AND  HA:  OON-HoRNE  KXtEFIMFNT  DAl'A 

No  two  experiments,  or  even  twe  sensors  fre.m  the  same  experiment, 
are  processed  in  exactly  the  .same  way.  Kowevei ,  the  p-r  rce.ssi  ng  of  any 
set  of  experiment  data  does  follow  a  series  of  common  steps.  This  sec¬ 
tion  describes  the  general  methodology  used  by  SA.SC  to  reduce  and  rrocess 
data,  regcirdless  of  the  specific  nature  of  the  sensor. 

Of  primary  importance  in  the  successiul  r'ro'cessinu  of  exper im.ent 
data  is  a  thorough  understanding  of  the  purtKise  of  the  experiment,  font 
of  the  data,  intermediate  stages  in  the  processing  stream,  anticipated 
results,  and  formats  for  data  presentation.  Discussion  between  SUWA  and 
SASC  of  the  processing  details  often  eliminates  delays  ana  makes  txissible 
a  flexible  system  design. 

Flexibility  is  of  extreme  im,oortance  because  the  processing  is  in 
many  cases  iterative,  with  modifications  at  each  .step-.  At  the  task  defi¬ 
nition  stage,  the  researcher  may  not  be  able  to  anticiiMte  the  full  exten.t 
of  the  processing  required.  Analysis  of  iiitermed i ate  viutputs  m.ay  bc- 
necessary  before  the  next  stage  of  protessiix;  cun  be  defined.  Cor.sc- 
quently,  algoritlims  must  be  developed  .in  ’,  ttie  coi  resixuviin  ;  software  de¬ 
signed  to  facilitate  implementation  of  the  lat  est  modi f  ications . 
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The  three  standard  pre-processing  information  inputs  are  the  analog 
strip  chart,  Digital  Work  Request,  and  Rocket  Data  Analysis  Information. 

The  analog  strip  chart  is  generated  at  the  rocket  launch  site  or  by  the 
Deconunutation  Section  of  the  AFGL  Computation  Branch.  The  strip  chart 
provides  the  experimenter  with  evidence  of  the  success  or  failure  of  the 
mission.  It  gives  the  data  processor  information  on  timing  or  noise 
problems  in  the  data.  Periods  of  extensive  noise,  which  requires  special 
filtering,  can  be  easily  identified  on  the  strip  chart. 

Figure  1  is  an  example  of  a  strip  chart  of  pre-flight  calibrations, 
showing  four  voltage  channels  (odd  numbers)  and  four  current  channels 
(even  numbers).  The  signal  is  in  telemetry  volts,  ranging  from  a  lower 
band  edge  of  zero  volts  (top  of  each  plot)  to  an  upper  band  edge  of  five 
volts.  In  the  calibration  stage  the  signal  is  stepped  0,  1,  2,  2.5,  3, 

4,  5  volts  and  then  back  to  zero.  Figure  2,  an  example  of  in-flight 
data,  provides  the  analyst  with,  the  opportunity  to  scan  the  transmitted 
data  for  any  instrument  problems. 

The  Digital  Work  Request  (Figure  3)  contains  information  about  the 
packed  data  that  are  to  be  processed.  For  example,  the  Figure  3  request 
form  shows  that  there  are  two  files  on  tape  R899  and  there  are  eight 
channels  (IRlG's  13  -  6).  The  short  time  on  file  1  (05:45:13  to  05:45: 

31)  indicates  that  it  is  a  calibration  file.  Other  information  contained 
in  the  Digital  Work  Request  includes  rocJcet  identification,  launch  date, 
sampling  rate,  start  and  stop  times  of  each  file,  and  data  frame  format. 

The  Rocket  Data  Analysis  Information  sheet  contains  a  description  of 
the  experiment  and  analysis  required,  calibrations  of  telemetry  volts  to 
engineering  units,  la'unch  date  and  time,  vehicle  number  and  type,  digital 
tape  format,  and  requested  items  such  as  listings,  tapes,  and  plots.  Fig¬ 
ure  4  is  an  example  of  page  one  of  the  Rocket  Data  Analysis  Information. 
Sometimes,  the  project  scientist  will  attach  technical  reports  describing 
the  instruments  involved  in  the  flight. 

Along  with  the  standard  pre-processing  inputs,  SUWA  provides  SASC  with 
the  appropriate  R-tape(s).  This  tape  contains  experiment  data  which  were 
telemetered  in  real  time  to  a  ground  recording  station.  A  typical  rocket 
flight  lasts  between  300  and  300  seconds;  the  sampling  rate  is  usually 
1,000  samples  per  second  (SPS)  but  may  go  to  over  10,000  SPS.  Data  on 
the  R-tape  are  in  packed  twelve-bit  words.  In  order  to  process 
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Digital  Work  Request. 


Problem  No.  3045 
Project  No.  76701010 

ROCKET  DATA  ANALYSIS  INFORMATION 

EXPERIMENTER:  Robert  R.  O’Neil _ 

EXPERIMENT  TITLE:  EXCEDE  -  SPECTRAL:  ELECTRON  ACCELERATOR 

LAUNCH  SITE:  POKER  FLAT,  ALASKA _ 

LAUNCH  DATE:  19  October  1979 _ 

File  1:  05:45:13  05:45:31 

LAUNCH  TIME:  05:46:40  AOS:  File  2:  05:47:40  LOS:  05:53:50 _ 

VEHICLE  NO:  A51.970 _  TYPE:  TALOS -CASTOR _ 

LINK  FREQ:  2289.5  MHz _ 

IRIG 

SUBCARRIER  ASSIGNMENTS:  13  ,  12  ,  11  ,  10  ,  9  ,  _8 _ ,  1,  6 

TYPE  OF  DATA;  PAM  _  FRAME  SIZE _ 

PCM  _  FRAME  SIZE  _ 

X  DIGITIZATION  RATE  1  HHz/  Channel 

PAM  OR  PCM  CAL  POINTS  _ ,  _ 

FM  INFLIGHT  CALS;  YES  _  NO  X 

PRELAUNCH  CAL.  DATA:  YES  X  NO  _  File  1  05:45:13  -  05:45:31 

ASPECT  REQUIREMENTS:  Yes  -  Separate  Request 

SPECIAL  REQUIREMENTS:  Accelerator  performance  to  be  merged  with  othtr 
EXCEDE  Measurements  (sensors),  aspect,  trajectory 

Start  Stop 

REQUIRED  START  TIME:  File  1  05:45:13  -  05:45:31 

STOP  TIME:  File  2  05:47:40  05:53:50 

Figure  4.  Rocket  Data  Analysis  Information. 
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these  data  on  the  CDC  6600  computer,  the  R-tape  is  unpacked  into  60-bit 
CDC-compatible  floating  point  words.  Unpacking  is  accomplished  by  using 
a  standardized  software  routine  designated  DQ  for  frequency  modulated  (FM)  data. 
When  the  data  are  pulse  code  modulated  (PCM)  a  standard  routine  HON316 
is  utilized. 

Once  the  data  are  unpacked,  quality  checks  are  performed  on  the  out¬ 
put.  SASC  uses  software  which  will  flag  any  two  consecutive  sampling 
times  which  are  separated  by  a  time  greater  than  the  sampling  rate,  or 
any  step  backwards  in  time.  Time  jumps  generally  come  in  pairs;  that  is, 
a  time  jump  backwards  is  usually  matched  in  magnitude  by  a  forward  time 
jump.  This  problem  results  from  the  software  used  to  generate  the  R- 
tape  and  can  usually  be  corrected  by  the  analyst.  However,  occasionally 
the  magnitude  of  the  backward  jump  does  not  equal  the  forward  jump.  It 
may  then  be  necessary  to  develop  an  algorithm  to  smooth  the  times  or  even  to 
create  a  new  R-tape.  When  the  timing  problems  have  been  resolved,  the  data 
are  put  on  tape  in  an  irreducible  form  and  the  specified  processing  begins. 

The  unpacked  data  tape,  without  any  timing  problems,  is  the  primary 
data  base.  The  information  on  this  tape  is  usually  in  counts  which  can 
later  be  translated  into  telemetry  volts  and  subsequently  into  engineering 
units.  Since  the  R-tape  is  returned  to  the  archives,  the  irreducible 
data  tape  serves  as  the  primary  source  of  experiment  data  and  remains  as  a 
point  of  origin  if  problems  are  encountered  downstream  in  the  data  reduc¬ 
tion. 

If  the  equation  which  converts  counts  to  telemetry  volts  is  known, 
a  new  data  base  may  be  created  immediately.  However,  if  the  calibration 
is  not  known,  it  is  necessary  to  list  data  in  the  vicinity  of  the  in¬ 
flight  calibrations.  Ihe  on-board  telemetry  system  is  able  to  send  a  step- 
pulsed  signal  (0  -  2.5  -  5.0  volts).  The  normal  lower  band  edge  is  repre¬ 
sented  by  a  count  value  of  -1638,  which  corresponds  to  zero  telemetry  volts, 
and  the  normal  upper  band  edge  has  a  count  of  1638,  which  translates  into 
five  volts.  The  analyst  averages  the  values  at  the  upper  and  lower  band 
edges  and  also  at  the  middle  of  the  band,  2.5  volts  (zero  counts).  Usually 
there  is  a  slight  shift  in  the  calibration  scale,  on  the  order  of  70  counts. 

A  linear  equation  is  derived  from  the  sample  listing  and  used  to  create  a 
data  base  in  telemetry  volts.  A  sample  listing  of  counts  and  corresponding 
telemetry  volts  is  made  and  an  analysis  is  done  to  check  the  validity  of 
the  derived  calibration  equation.  In  most  cases  a  fine-tuning  of  the  cali- 
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bration  equation  will  yield  the  appropriate  telemetry  voltages.^ 

When  the  Contract  Manaaer  is  confident  that  the  appropriate  linear 
equation  to  convert  counts  into  telemetry  volts  has  been  found,  the  data 
are  ready  to  be  changed  into  engineering  units.  The  calibration  equations 
used  to  give  physical  meaning  to  the  count  data  received  at  the  ground  sta¬ 
tion  are  the  result  of  work  by  the  instrument  technicians  in  their  labora¬ 
tory.  The  simplest  eouation  is  linear  and  valid  for  any  value  of  the  tele¬ 
metry  voltage.  The  more  complicated  equations  involve  logarithmic  scales 
and  may  be  valid  only  over  a  st^ecified  range. 

Computer  application  of  the  calibration  equation (s)  to  the  raw  data 
produces  a  data  base  on  magnetic  cape  with  information  in  units  that  have 
physical  meaning.  This  data  base  is  ready  to  be  reduced,  listed,  or  dis¬ 
played.  If  the  sampling  rate  is  1,000  SPS ,  it  is  usually  sufficient  to 
list  the  data  at  20  -  50  SPS.  This  avoids  unnecessary  and  extensive  list¬ 
ings  which  take  up  computer  time  and  space.  If  the  experiment  dictates  a 
high  Scimplinq  rate  at  specified  times,  it  becomes  a  matter  of  retrieving 
the  data  from  the  original  engineering  unit  tape. 

Data  may  be  displayed  by  pen-and-ink  plots,  microfiche,  or  Tektronix. 
Pen-and-ink  plots  are  best  suited  for  later  reproduction  in  technical  pub¬ 
lications.  Microfiche-  plots  permit  the  most  plots  in  the  least  amount  of 
space.  Copies  of  microfiche  plots  are  suitable  for  publication  purposes 
but  yield  to  pen-and-ink  because  of  .slight  distortions  in  the  reproducing 
process.  Plots  which  are  copied  directly  from  the  Tektronix  are  useful  in 
the  preliminary  stages  of  p'lotting  the  data,  because  the  turn-around  time 
is  short. 

>^103117,  after  the  data  have  been  reduced,  trajectory  and  attitude 
data  are  merged  with  the  experiment  data.  Standard  trajectory  and  attitude 
tapes  are  provided  to  SASC  to  create  a  final  tape  containing  all  the  infor¬ 
mation.  If  listings  or  plots  of  the  combined  data  are  required,  they  are 
supplied. 

C.  DATA  PRCX;ESSING:  FLIGHT  B  SPECTROMETER 

Although  there  is  a  wide  range  of  experiments  at  AFGL,  the  proces¬ 
sing  of  rocket-  and  balloon-borne  sensor  data  follows  a  series  of  common 

Vor  cl  detailed  study  of  telemetry  systems  the  reader  is  encouraged  to  re¬ 
fer  t^'  Aerospace  Telemetry,  Vol.  1,  edited  by  Harry  L.  Stiltz.  Prentice- 
Hall  Inc.,  Englewood  Cliff.s,  New  Jersey. 
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steps.  This  section  describes  the  step-by-step  process  followed  in  an 
actual  project. 

Scientists  at  AFGL  are  conducting  a  series  of  experiments  to  obtain 
absolute  values  of  the  solar  flux  over  a  rising  period  of  the  eleven-year 
solar  cycle.  Measurements  are  made  in  the  vacuum  and  near  ultraviolet 
regions.  Five  rocket  flights  were  made  from  May  1976  to  September  1980, 
as  listed  in  Table  1. 


TABLE  1.  FIVE  ROCKET  LAUNCHES  -  SPECTROMETER 


FLIGHT  DATE 


VEHICLE  AND  INSTRUMENTATION 


A  May  18,  1976 


B  April  21,  1977 


C  August  9,  1977 


D  September  19,  1978 


E  September  18,  1980 


Rocket  -  A03. 410-1 

Instrument  -  Rocket  Spectrometer  No.  61-1 
Wavelength  range  -  170  -  350  nm 
Number  of  Spectrometers  -  1 

Rocket  -  A04.410-2 

Instrument  -  Rocket  Spectrometer  No.  61-II 
Wavelength  range  -  170  -  350  nm 
Number  of  spectromters  -  1 

Rocket  -  A03. 509.1 

Instrument  -  Rocket  Spectrometer  No.  62-1 
Wavelength  range  -  130  -  180  nm;  170  - 
350  nm 

Number  of  spectrometers  -  2 
Rocket  -  A04.711-1 

Instrument  -  Rocket  Spectrometer  No.  62-11 
Wavelength  range  -  130  -  180  nm;  170  - 
350  nm 

Number  of  spectrometers  -  2 
Rocket  -  A04.801(RS-61) 

Instrument  -  Rocket  Spectrometer  No.  62-III 
Wavelength  range  -  175  -  310  nm 
Number  of  spectrometers  -  1 


This  section  will  cover  Flight  B. 


At  the  initial  meeting  with  the  SUWA  Contract  Manager,  SASC  was  given 
a  copy  of  the  Digital  Work  Recfuest  (Figure  5) ,  which  described  the  R-tape 
R730,  and  a  detailed  memorandum  which  outlined  the  data  processing  require¬ 
ments. 
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The  raw  telemetry  data  for  each  fliqht  were  generated  by  a  solar- 
pointed  spectrometer.  The  spectrometer  admits  light  through  a  narrow 
entrance  slit,  separates  the  light  according  to  wavelength  by  an  optical 
grating,  and  focusses  the  image  of  the  slit  on  the  photosensitive  surface 
of  a  detector  which  in  turn  produces  an  electrical  signal  proportional  to 
the  incident  illumination.  This  detector  signal,  called  photon  count,  is 
transmitted  digitally  in  4  BCD  bits  which  range  from  0000  to  9999. 

(On  a  logarithmic  plot  0  is  to  be  plotted  as  1,  (10*^)  ,  and  hence  a  four- 

decade  axis  is  required  to  cover  the  full  range  of  data,  namely,  10*^, 

12  3 

10  ,  10  ,  10  .) 

Wavelengths  are  selected  by  mechanical  rotation  of  the  optical  grat¬ 
ing-  This  rotation  is  transmitted  through  a  gear  train  from  a  pulsed 
stepping  motor,  connected  to  an  encoder.  The  encoder  provides  a  4-bit 
BCD  number  to  telemetry  to  establish  the  wavelength  at  which  the  photon 
count  data  were  recorded.  The  encoder  does  not  have  sufficient  resolu¬ 
tion  to  indicate  the  motion  associated  with  one  drive  pulse,  but  rather 
increments  once  every  13  steps.  The  telemetry  data  will  therefore  show 
the  same  wavelength  value  13  times  for  a  sequence  of  data  frames,  although 
the  wavelength  has  in  fact  changed.  The  standard  time  interval  between 
frames  is  4  ±  1  milliseconds. 

The  first  step  in  the  data  analysis  and  reduction  was  to  unpack  the 
twelve-bit  words  into  60-bit  CDC  words.  Tape  R730  was  unpacked  onto  tape 
CC0487  and  returned  to  the  AFGL  tape  archives.  Tape  CC0487  became  the 
primary  data  source,  arranged  in  digitized  format.  Figure  6  shows  the  bit 
construction  of  the  words  on  CC0487.  Each  frame  of  data  contains  four 
words,  including  universal  time  of  day  in  seconds,  wavelength,  and  pho¬ 
ton  count.  As  an  example,  consider  a  word  which  has  a  zero  fill  for 
bits  60  to  17,  and  bits  16  to  1  look  like: 

Bit  16  15  14  13  12  11  10  9  8  7  6  5  4  3  2  1 
0  1  0  0  0  0  010011010  1. 

This  word,  when  divided  into  four  different  values,  each  multiplied  by  an 
appropriate  power  of  ten ,  yields ; 

5(10°)  +  3(10^)  +  I(IO^)  +  4(10^)  =  4135. 

The  next  step  was  to  create  a  data  base  which  allowed  the  analyst 
wide  latitude  in  processing  and  reduced  computer  time.  This  new  data 
base  was  created  from  CC0487  and  put  on  tape  CC3901.  Each  data  frame  con- 
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60  bit 

photon  count  word 


60  bit 

time, of  day 
(sec) 


60  bit 

60  bit 

sync  word 

wavelength  word 

Each  has  16  bits,  right  justified,  with  a 
zero  fill 


WAVELENGTH  WORD  (bit  60  to  17  zero-filled) 


xlO  xlO  xlO  xlO 

Bit^6~i5*i4  13^l2~il  ^0  9  87654321 


PHOTON  COUNT  WORD  (bit  60  to  17  zero-filled) 


Bit  16  15  14  13  12  11  10  9  8  7  6  5  4  3  2  1 


Figure  6. 


Bit  construction  of  words. 


sisted  of  universal  time  of  day  (HHMMSS) ,  integer  wavelength  (angstroms) , 
and  integer  photon  count.  A  listing  of  these  raw  data  was  also  created 
and  delivered  to  the  Contract  Manager.  Figure  7  is  a  sample  of  that  list¬ 
ing  . 

As  noted  above,  the  mechanical  arrangement  of  the  instrument  dictates 
13  steps  per  wavelength  and  sampling  at  4  ±  1  milliseconds.  After  tape 
CC3901  was  created  it  was  necessary  to  develop  an  algorithm  which  would 
flag  all  deviations  from  the  instrumentation  routine.  The  Contract  Mana¬ 
ger  provided  SASC  with  the  time  segments  involved  in  the  spectrometer  scans, 
excluding  the  time  when  the  filter  is  inserted.  The  data  on  CC3901  were 
checked  for  irregular  wavelength  and  time  intervals.  All  irregularities 
were  identified  by  SASC  and  this  information  was  supplied  to  the  Contract 
Manager.  Recommendations  for  adjusting  the  data  base  to  correct  for  the 
irregularities  were  then  given  to  SASC. 

The  changes  to  CC3901  were  made  and  a  complete  data  base  with  regu¬ 
lar  wavelength  and  time  sampling  was  made  on  CC4521.  A  quality  check  of 
data  changes  was  made  and  the  entire  data  base  was  listed,  in  the  format 
of  Figure  7,  and  displayed  on  microfiche,  as  in  Figure  8.  At  this  point 
the  data  were  correct,  according  to  the  sampling  rate,  but  errors  in  the 
count  rate  caused  by  bit  transmission  difficulties  had  to  be  corrected. 

Identification  and  correction  of  the  transmission  errors  can  be 
done  by  hand  but  is  a  time  consuming  process.  An  algorithm  was  developed 
as  a  result  of  several  meetings  of  the  SASC  Principal  Investigator,  Con¬ 
tract  Manager,  and  project  representative.  The  full  algorithm  (see  Section 
C.l)  was  carried  out  on  the  CDC  6600.  The  ensuing  refined  data  base  was 
listed  with  the  bad  data  (in  parentheses)  replaced  with  corrected  data 
(Figure  9)  and  put  on  tape  CC1520. 

At  this  point,  80  -  90%  of  the  transmission  errors  had  been  identi¬ 
fied  and  adjusted.  The  remaining  10  -  20%  were  identified  by  the  project 
team  and  listed  for  SASC.  Software  was  developed  to  modify  the  data  base 
on  CC1520.  A  final  data  base  free  of  both  sampling  and  transmission 
errors  was  made  on  CC1974,  and  a  corresponding  data  list  was  given  to  the 
Contract  Manager. 

The  next  step  will  be  to  display  the  data  on  pen-and-ink  plots.  When 
the  project  is  completed,  SASC  will  deliver  a  final  data  tape  and  all  soft¬ 
ware  developed  by  SASC  which  is  determined  to  be  useful  in  the  future. 
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Figure  8.  Sample  of  microfiche  plot  of  Flight  B  data. 
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RECO  FLIGHT  B  APRIL  21»1977 
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1. 


ERROR  REMOVAL  METHOD 


Error  Identification  in  the  Data  Set  n, ,n_,....n. 

12  N 


Definitions: 

(a)  n  ,  n,  ,  n  ,  n^  -  four  consecutive  values  in  the  data  set. 
abed 


-  "b 


An,  =  n,  -  n 
be  D  c 

An  ,  =  n  -  n, 
cd  c  d 


(c)  n 


ab 


n 

a 

+ 

"b 

T 

"b 

+ 

n 

c 

T 

n 

c 

+ 

"d 

cd 


(d)  G 


(e)  PCK 


An  .  An. 
ab  ,  be 

S  ■  "ab  "be 

An^  An  , 
^  _  be  cd 

=  — = —  •  -= - 

C  n,  n  , 

be  cd 

1 


ac 


Iterations  : 


(a)  The  first  pass  through  the  data  set  makes  changes  to  errors 
which  are  the  result  of  transmission  problems  with  the 
1,000  word.  An  error  is  obtained  when 


ab 


n,_  <  -  10 

be 


(b)  The  second  pass  through  will  correct  errors  which  may  range 
from  2  to  880. 


Identification  of  errors. 

(a)  First  pass:  n  .  •  n,_  <10 

ab  be 

(b)  Second  pass: 
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(i)  Compute  G  .  If  G  i  0  there  is  no  error,  but 

B  B  - 

if  G  <  0  proceed  to  (ii) . 

B 

(ii)  Compute  PCK.  iflc  |  <  PCK  there  is  no  error; 
otherwise  proceed  to  (iii). 

(iii)  Compute  |g^|*  If  >  jG^jthe  error  is  at  b; 

otherwise,  error  is  at  c. 

4.  Correction  of  error.  Consider  three  consecutive  data  points 
x^.x^,  x^  where  the  error  has  been  identified  at  x^. 

(x^  +  x^) 

(a)  Compute  average.  AV'E  =  - 2 - 

(b)  Compute  the  absolute  difference  between  the  error  and  AVE. 

GLDIF  =  jxj  -  AVe] 

(c)  Calculate  the  difference  between  GLDIF  and  each  factor. 

The  first  pass  uses 

FACTOR  (j)  j  =  1,..,4  which  are 

1,000,  2,000,  4,000,  8,000  respectively. 

A  list  of  the  38  factors  used  for  the  second  pass  is  in¬ 
cluded  at  the  end  of  this  section. 

DIF  (j)  =  [(GLDIF  -  FACTOR  (j))) 

(d)  Determine  which  FACTOR  (j)  is  the  closest  to  zero. 

(e)  Determine  whether  FACTOR  (j )  is  added  to  or  subtracted 
from  the  error  x^. 

5.  Listing  of  data. 

(a)  Errors  -  The  time  and  wavelength  associated  with  the 
error,  the  error  and  its  new  value,  and  the  measurement 
before  and  after  the  error  are  listed. 

(b)  Data  -  The  final  listing  of  the  data  shows  time  -  wavelength 
count.  In  the  case  where  an  error  was  discovered,  the  pre¬ 
vious  value  of  the  count  is  provided  in  parentheses. 

(c)  Magnetic  tape  -  The  format  of  the  words  on  magnetic  tape 
is: 

IHR,  MIN,  SEC,  IWAVE,  KOUNTN,  KOUNTO,  IFLG 
where  IHR,  MIN,  SEC  are  hour s-minutes- seconds  (of  day) 
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IWAVE  -  wavelength 

KOUNTN*  -  present  value  of  count 

KOUNTO  -  old  value  of  count 

IFLG  -  0  indicates  data  frame  is  unaltered 
1  indicates  data  frame  is  changed. 


KOUNTN  =  KOUNTO  in  those  frames  where  no  error  is  found 


6.  List  of  factors  for  second  pass: 


1. 

2 

2. 

4 

3. 

8 

4. 

10 

5. 

12 

6. 

14 

7. 

18 

8. 

20 

9. 

22 

10. 

24 

11. 

28 

12. 

40 

13. 

42 

14. 

44 

15. 

48 

16. 

82 

17. 

84 

18. 

86 

19. 

88 

20. 

20 

21. 

40 

22. 

80 

23- 

100 

24. 

120 

25. 

140 

26. 

180 

27. 

200 

28. 

220 

29. 

240 

30. 

280 

31. 

400 

32. 

420 

33. 

440 

34. 

480 

35. 

820 

36. 

840 

37. 

860 

38. 

880 
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